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Abstract 

Obtaining ultracold samples of dipolar molecules is a current challenge which requires an ac- 
curate knowledge of their electronic properties to guide the ongoing experiments. In this paper, 
we systematically investigate the ground state and the lowest triplet state of mixed alkali dimers 
(involving Li, Na, K, Rb, Cs) using a standard quantum chemistry approach based on pseudopo- 
tentials for atomic core representation, gaussian basis sets, and effective terms for core polarization 
effects. We emphasize on the convergence of the results for permanent dipole moments regarding 
the size of the gaussian basis set, and we discuss their predicted accuracy by comparing to other 
theoretical calculations or available experimental values. We also revisit the difficulty to compare 
computed potential curves among published papers, due to the differences in the modelization of 
core-core interaction. 
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I. INTRODUCTION 



Researches on ultracold molecules are progressing very fast since a couple of years, as a 
growing number of research groups are now involved in this field. Indeed, ultracold molecules 
offer many new opportunities compared to ultracold atoms, due to their more complex struc- 
ture: they can carry an amount of internal energy larger by orders of magnitude than their 
rinetic energy, opening new areas like ultracold photochemistry |l| or also superchemistry 
3, The ultimate control of elementary chemical reactions using short-pulse lasers is now 
foreseen as both the internal and external degrees of freedom of the molecules can be 
well mastered. They allow the reinvestigation of phase transitions which up to now were 
explored only from the point of view of condensed matter physics, like the transition from 
bosonic to fermionic statistics , the superconductivity transition through the obser- 

vation of the pairing gap in strongly interacting Fermi gases |§] or the superfluid transition 
via the study of the BEC-BCS crossover Spectacular achievements have been recently 
demonstrated as for instance the observation of molecular condensates with various alkali 
dimers 10 , 11 , 12 ] , or of Feshbach resonance in molecule-molecule ultracold collisions Q| . 

Other interesting prospects concern dipolar molecules (ie with a permanent dipole mo- 
ment), as it can be read through the special issue on this topic published in 2004 [rsj ]. 



Among those are the possibility to design a quantum information device |l£ 



or the test of 



fundamental theories like the measurement of the eletron dipole moment |l6l. Il7|. 

The formation process is obviously the key issue of these developments. An efficient 
approach which dramatically improved these last years is the so-called Stark deceleration 
technique, relying on the slowing and trapping of molecules with a permanent electric dipole 
moment (usually labelled as dipolar molecules) using inhomogeneous external electric fields 
Q| . Besides, the " 'historical" ' approach based on photoassociation of ultracold atoms which 
has first demonstrated the formation of ultracold molecules with cesium dimers |19], has 
recently progressed with the creation of ultracold samples of dipolar molecules composed of 



different alkali atoms like RbCs |2ij,|2l|, KRb |22j,|2j, and NaCs |2J,|25(. 

In a previous paper we computed the rates for the photoassociation of mixed alkali 
pairs, and for the susbsequent formation of cold molecules, which showed that all alkali pairs 
involving either Rb and Cs are well suited for that purpose, as the cold molecule formation 
rate was only about ten times smaller than for Cs2 formation. This result was found in good 
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agreement with the RbCs experiment of Kerman et al |20|. Potential curves available in the 
literature, and constant (atomic) transition dipole moment were proved to be sufficient to 
establish these estimates. However, the practical implementation of cold molecule formation 
via photoassociation will require a much better knowledge of their electronic properties 
like radial variation of permanent and transition dipole moments, in order to better guide 
the experimentalists towards specific systems and transitions with maximal efficiency. New 
spectroscopic studies on heteronuclear alkali pairs are also currently reinvestigated for NaRb 
27 1 and NaCs 2^|, as well as the investigation of the effect of strong electric fields 2^1, which 
may also benefit of such a study. 

To respond to these requirements, we started a new accurate analysis of electronic prop- 
erties of all alkali pairs from Li to Cs, including potential curves for ground and excited 
states, permanent and transition dipole moments. We set up an automatic procedure based 
on the CIPSI package (Configuration Interaction by Perturbation of a multiconfiguration 
wave function Selected Iteratively) [3| developed by the "Laboratoire de Physique Quan- 
tique de Toulouse (France)". Here we first investigate the permanent dipole moment of 
the ground state and the lowest triplet state of all mixed alkali pairs, for which only scat- 
tered experimental or theoretical results are available. We display their variation with the 
interatomic distance, as well as with the vibrational level. Most of these results were not 
previously available elsewhere. We emphasize on their convergence with respect to the size 
of the basis set. In the next section, we briefly recall the calculation procedure, which is well 
documented in many publications. We also review extensively the available calculations on 
these systems. In section UTTj we present our results obtained with basis sets published in the 
literature, and with new extended basis sets designed for all alkali atoms, in order to address 
convergence issues with respect to the size of the basis. We also discuss the comparison with 
other published papers, which is not often proposed, especially in the perspective of the 
addition of effective terms to the computed energies. 



II. THEORETICAL CALCULATIONS 



Many studies on heteronuclear alkali dimers have been performed using various quantum 
chemistry methods, restricted to the calculation of electronic potential curves for the ground 
and excited molecular states. The CIPSI package has been used to compute ground and 
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excited potential curves of RbCs 



NaCs, KCs |38 
and Barthelat 



31 



3IQ, NaK 3,0, KLi 3> Li 



LiRb, NaRb 



3J, LiCs, 



KRb 39]. The method is based on the ab initio pseudopotentials of Durand 



40l . |4l| for core representation, phenomenologic ^-dependent effective core 
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polarization potential |42| and self-consistent calculation (SCF) combined with full valence 
configuration interaction interaction (CI) calculations. Additional effective terms may also 
be used to take in account short range core polarization terms (see for instance ref. 
Potential energy curves of ground and excited states of LiNa have been investigated by 
Davies et al |43| through many-body perturbation theory calculations, and by Schmidt- 
Mink et al k4| who have combined all-electron (SCF) calculations, valence CI calculations 

n 

and have included an effective core polarization potential |45(. The electronic structure of 
KRb has been studied by Leininger et al and by Park et al 3] using a method based 
on restricted Hartree-Fock (RHF) calculations, using the small-core (with nine electrons for 
each atom) relativistic pseudopotential of Christansen et al j3| and CI calculations. Park et 
al also introduced the effective core polarization potential of ref. j^lj . Note that an extensive 
bibliographic overview on theoretical as well as experimental papers can be found in the 
"DiRef database 0] devoted to diatomic molecules. 

Only few papers include theoretical values for permanent dipole moments in addition 
to potential curve calculations. Igel-Mann et al |52| have performed a systematic study of 
the permanent dipole moments around the equilibrium distance of the ground state of all 



alkali dimers using the semi-empirical pseudopotential of Fuentealba |5J| including a core 
polarization potential. All-electron SCF and valence CI calculations with or without core 



polarization correction have been performed by Miiller and Meyer [45j] for LiNa, LiK and 
NaK at the equilibrium distance. Janoschek et al j^4| and Stevens et al 55] have reported 
permanent dipoles of NaK. The ab initio and pseudopotential calculations of Janoschek 
concern only the ground state, while the more complete SCF and CI calculations including 
core-polarization effects done by Stevens et al [5_5j include the .R-dependence of the ground 
and triplet state permanent dipole moment. SCF and CI calculations of permanent dipole 
moments for NaLi have been performed by Bertoncini et al and by Rosmus and Meyer 



571 ] . Bertoncini et al reported the -R-dependence of the dipole in the ground and the 



lowest triplet states. 

Ab initio calculations of the independent dipole moments of the 1,3 S + states of KRb have 
been performed by Kotochigova et al [58], using Hartree-Fock or Dirac-Fock orbital basis 
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set complemented by Sturmian functions to describe virtual orbitals and CI valence-bond 
calculations. Core-valence correlations are accounted by allowing single electron excitation 
from the closed 3p 6 of K and 4p 6 of Rb and by introducing several virtual orbitals. Core 
polarization effects are not explicitely mentionned by the authors in their papers. They 
found a large influence of relativistic effects on the dipole moments, with an increase of 
about 20% for the ground state and of about 50% for the a 3 S + case, in the region of the 
potential minimum. 

Besides, only a few experimental values for the permanent dipole moments of the 
ground state obtained by different techniques have been reported. Dagdigian and Wharton 
59, 0, 1(3]] used molecular beam electric deflection and resonance spectroscopy to deter- 
mine the molecular electric dipole of several heteronuclear alkali dimers. The permanent 
dipole of NaLi or KLi molecules have been determined by Engelke et al using laser induced 



fluorescence I62l 1631 1 . and the one of NaK by Wormsbecher et al using microwave optical 

n n 

double resonance technique |64Jj . Finally, Tarnovsky et al |£5J] have measured the electric 
dipole polarizabilities of alkali homonuclear and NaK and KCs heteronuclear dimers using 
molecular beam deflection in an inhomogenous electric field. This led the authors to set up 
two empirical rules for determining polarizabilities and dipole moments of all heteronuclear 
alkali dimers, yielding then the values of the dipole moments for all mixed pairs from the 
difference between the polarizabilities of their respective homonuclear dimers. 

As mentionned above, we use here the CIPSI package to investigate the electronic struc- 
ture of mixed alkali pairs MM', and we recall now the main features of this code. In this 



met 
[40. 



rod, the atomic cores are described by the pseudopotentials of Durand and Barthelat 



41 1 whose parameters have been adjusted to reproduce the energies and valence orbitals 



of all-electron Hartree-Fock SCF calculation for the atomic ground state. For Rb and Cs, 
the pseudopotential includes the mass- velocity, and the Darwin relativistic corrections. Each 
atom is described by a Gaussian basis set. As the alkali atoms are treated as one-electron 
species, correlations are not explicitly introduced. Core polarization is taken in account 
through an ^-dependent effective core polarization potential (ECP) depending on the 
dipole polarizability a^ I+ and a^' of the M + and M' + ions and of I- dependent cut-off 



66 
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681 of the two lowest 



parameters pi adjusted to reproduce the experimental energies 
s, p and d atomic levels. Molecular orbitals are then determined by restricted Hartree-Fock 
calculations and full valence CI are performed. 
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Additionnally, the potential energy should include the core-core interaction potential 
V CC (R) which is modelled as long as the cores are far enough from each other with the 
pure repulsive 1/R term, and an attractive charge- induced dipole term V™ d (R) = — (ar^ + + 
a% r ' + )/2R A . When the cores come close to each other, such an approximation is not sufficient 
due to their strong electrostatic repulsion, and further short-range effective terms have to 
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69, 701, the 



be included in V cc . As discussed for instance in the case of alkali dimers 
core-core repulsion V^ P (R) has to be evaluated from ab-initio calculations according to 
various assumptions. The frozen core approach is then often used, in which the repulsion 
between the two ionic cores is obtained from a SCF calculation where the two ionic cores are 
represented by frozen atomic orbitals. It has been shown jy| that this term can be fitted 
by an exponential form. The contribution of the core-core dispersion energy has also to be 
taken in account, approximated by the London formula (see for instance ref. 

ydis P(m _ 6a d a d h i h i m 
cc K ' 2RS Ef I+ +Ef I,+ { ) 

where Ef + and Ef + are the ionisation energies of the M + and M' + ions respectively. 

We performed three different series of calculations, described below, characterized by 
different Gaussian basis sets summarized in Tables |Hj and IIHI The resulting atomic 
energies are displayed in Table IIVI while the size of the generated molecular basis set is 
given in Table IVI for mixed alkali pairs. 

• Prior to our new calculations, we checked the code and our procedure by first repro- 
ducing the computations performed by previous authors, employing the contracted 
Gaussian basis sets, the ionic core dipole polarizabilities and the ^-dependent cut-off 
radii pi published by Poteau et al jzij for Li, by Magnier et al 7^| for Na, by Magnier 
et al 34] for K, and by Pavolini et al for Rb and Cs 1311. Note that in these calcu- 

nn 

lations, the values of may come from different sources |45|, 73]. These calculations 
will be refered to as belonging to the "A" series. 

• Then we removed the contracted orbitals from the basis sets quoted just above, and 
introduced the ^-dependent ECP's afterwards. For all atoms, we obtained a better 
description for atomic energy levels than with the contracted basis sets. Guided by 
this trend, our next series of calculations (labelled as series " B" ) involve uncontracted 
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Gaussian basis extended with respect those used in the "A" series, and dipole polariz- 
abilities all coming from the same paper by Wilson et al Some basis coefficients, 
and the cut-off radii, have been reoptimized, in order to improve the quality of calcu- 
lated atomic energy levels. Note that the Li basis set has been considerably increased 
compared to ref . |7ll| . In contrast, the K basis set is the same than in ref.j3|, but 
without contraction coefficients. The K basis has been recently extended by Magnier 
et al [l^ to study highly excited states of the K 2 dimer. 

Due to the current interest for cold molecule studies about heteronuclear alkali pairs 



involving Cs, 0, 



25I ]. we set up a third series of computations (C series) for LiCs, 



NaCs, KCs, and RbCs, using a basis set and cut-off parameters for Cs, which will be 
illustrative for the discussion of convergence issues. 

As seen in Table HV1 we obtained in most cases slight improvements on the atomic energy 
levels when A and B calculations are compared to the experimental ones, which confirms 
that previously published basis sets were already quite well optimized. However, as it will 
be discussed in the next section, the improvement is more spectacular for potential curves, 
due to the cumulative effect of better defined static polarizabilities, and extended basis sets. 



III. RESULTS 



Before discussing in detail our results on permanent dipole moments, we first check the 
reliability of our extended basis sets by looking at the potential curves obtained for the 
ground state and lowest triplet state of the mixed pairs, concentrating on two specific 
examples, RbCs and KRb. We first first look at the potential energies V'(R) including 
only the 1/R term of V CC (R). As it is shown in Figure ^ the minimum energy for V'(R) 
curves is found in general at distance R m larger than the equilibrium distance R e of the full 
V(R) = V'(R) + V™ d curves. The combined influence of the addition of diffuse orbitals and 
the improved adjustement of the cut-off radii in series B compared to series A is clear, as 
the X and a states are found 48 cm -1 and 6 cm" 1 deeper than in series A, respectively. In 
the latter case, the basis C for Cs still increases the a well depth by 6 cm -1 , while no effect 
is visible for the X state. The same trend is observed for all pairs: as a further example, the 
KRb ground state is found in the B case about 70 cm -1 deeper than in the A case, while 
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the depth of the a state is increased by 2 cm -1 . Excited potential curves are also modified 
from A to B calculations, and will be investigated in further work. 

We then computed the permanent dipole moment of X and a states for all heteronuclear 
pairs. The sign of the permanent dipole moment of a molecule depends on the interatomic 
axis orientation, which is assumed here to be oriented along the MM' direction, where M 
is the lightest atom of the mixed pair. A negative dipole moment then implies an excess 
electron charge on the M atom. In the following, values are given in Debye, with 1 atomic 
unit=2.54158059 Debye. 

Here again we first check the influence of the optimization of the basis set and of the 
cut-off parameters on the RbCs permanent dipole moment (figure |2J). Note that V CC (R) does 
not influence the electronic wave functions, neither the dipole moments functions. The R 
dependence is found very similar for all A, B, C calculations, with a minimum around 9.3ao 
for the X state (slightly larger than the minimum distance R^ = 8.5%), and around 13.4a 
for the a state (slightly larger than the minimum distance R^ = 12.6cto). Increasing the size 
of the basis changes the dipole moment of the X state in the region of the minimum by only 
2 to 5%, while the moment for the a state, which is much smaller in magnitude, changes 
(around the location of its minimum) by about 50% from A to B calculations, and only by 
7% from B to C calculations. A similar trend is observed for the other mixed pairs (see 
table IVl]) . as the change from A to B or C approaches in the ground state dipole moment 
never exceeds a few percents, being most often smaller than 1%. So we can safely estimate 
that the B calculations are converged with respect to the size of the basis set. 

The R- variation of the permanent dipole moment for the A 1 S + and a 3 E + states deduced 
from calculations B (figure EK) is found similar for the ground state of all pairs. Apart for 
the LiNa dipole moment which remainsclose to Debye, the curves present a minimum at 
a distance systematically larger by about la to 1.5ao than the equilibrium distance 
Rf of the ground state, and vanishes at large R as expected. The magnitude of the dipole 
strongly depends of the chosen pair, the largest being for LiCs around -6 Debye. In contrast, 
the independence of the dipole moment in the a state is more irregular over the mixed pair 
series, with a minimum for some species (LiCs, RbCs, KRb) and a maximum for the other 
dimers. 

Figure 0] shows the variation of the permanent dipole moment with the vibrational level of 
the X and a states, obtained after averaging the dipole moment function over the vibrational 
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wave functions of the potential functions V(R) = V'(R) + V™ d . This figure partly reflects 
the .R-dependence drawn in Figure EJ the ground state dipole moment is almost constant for 
the lowest forty levels for the less polar species, while it varies more rapidly for LiRb, LiCs, 
NaRb. As expected, all these curves vanish at large distances (not shown on the figure). 

For the purpose of comparison with other works, we summarize our results for the ground 
state in a tabular form (Table ETJ), displaying the value of the permanent dipole moment 
at the equilibrium distance R* , at its extremum value R^, and averaged for the v — 
vibrational level. To the best of our knowledge, the other available theoretical values for 
the permanent dipole moment of the mixed alkali pairs concern the ground state at its 
equlibrium distance Rf . In contrast, only the value of the permanent dipole moment for 
a given vibrational level of molecules is obtainable in an experiment, limited to v — 
up to now. Note that for simplicity we reported absolute values in table IVII indeed, the 
experimentalists cannot determine the sign of the dipole moment, while it is not always 
clear in the theoretical papers what is the choice of the authors for the orientation of the 
molecular axis. 

As expected, the dipole moment of v = is close to the value of the dipole function at 
the equilibrium distance. Our results are in good agreement with the available experimental 
values, within a 2% for LiK, NaK, and NaCs, and slightly more (6%) for NaRb. The 
difference is found much larger (around 20%) for the less polar molecule LiNa. We confirm 
that the approximate scaling law proposed in ref. overestimates the dipole moment, 
except for the less polar molecules KRb and RbCs, which in contrast are underestimated. 

n 

Our results are in excellent agreement with those of ref. [52] , which a similar approach without 
the ^-dependence of the core polarization terms. 

We have already mentionned that very few other theoretical variations of the perma- 
nent dipole moment with the interatomic distance are available in the literature. The 
-R-dependence of the dipole moment in the ground state of NaLi molecule obtained by 

n n 

Bertoncini et al [56] and by Rosmus and Meyer |5JZJ disagree with each other both in sign 
and in magnitude. Our negative value for LiNa supports the prediction of Rosmus and Meyer 
but not their R- variation. Note that the authors of ref.js2|, and later on in another paper 
5oj | , underlined the great sensivity of the dipole moment of NaLi to valence correlations and 
core polarization effects. Similarly, results by Stevens et al on NaK 55] also differ strongly 
from ours, which can be attributed to the absence of an ^-dependence in the effective core 
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potentials, included afterwards in the NaK study of Magnier and Millie 

The KRb molecule is an attractive case in the cold molecule field, as two recent calcu- 
lations of the X and a states of KRb includin g th e .R-dependence of the permanent dipole 
moments have been performed by Park et al |47. |3| and by Kotochigova et al 5^1 using 
quite different theoretical approaches. Their predictions are compared with our calculations 
A and B in figure El It is striking to see that quite large differences are observed, up to 20% 



with the non-re 
culations of ref. 



ativistic result of ref. 



581 ] . and 40% with ref. |47( or with the relativistic cal- 
581. It is not easy to understand such large differences. One reason could be 

n 

the different ways to treat the core polarization effects. Kotochigova et al |5J| indicate that 
the implementation of core polarization effects may reduce their dipole moment by about 
50%. They also mention that their convergence with respect to the number of basis func- 
tions is of the same order than the difference between their relativistic and non-relativistic 
calculations. It would be interesting to perform these kind of comparisons for molecules 
with a larger dipole moment, as we may expect larger core polarization effects. Given the 
good agreement we obtain with available experimental values on other systems, we believe 
that our values are quite accurate. 

IV. DISCUSSION 

We performed potential curves and transition dipole moments calculations for all het- 
eronuclear alkali pairs, in order to investigate the influence of the basis set size on higher 
excited electronic states, which is expected to be significant. The comparison of these poten- 
tial curves with other existing calculations will be presented in further publications, as well 
as results for dipole transition moments which are not available elsewhere for most species. 
These data are of particular importance for the optimization of multiple-step schemes for 
ultracold molecule formation, like the one performed by Sage et al j^jj, or for the interpre- 
tation of recent observation of Autler-Townes effect in highly excited molecular states 79] . 
An important issue is also the evaluation of the accuracy of potential curves (by looking 
for instance at well depths and equilibrium distances) compared to experimental determi- 
nations, in relation with the choice of the effective terms included in the core-core potential: 
Table IVTU displays the contributions of the V^ sp term, as evaluated with equation [Tl and of 

nn 

the V™ v term from refs. [3JJ, Hii| close to the equilibrium distance of the ground state of the 
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KRb and RbCs systems. The largest effect is obtained for the heaviest pair RbCs. Further- 
more, the V^ sp term should be cut off at a somewhat arbitrary distance, when ionic cores 

n 

become close to each other. Jeung |69| discussed various approximations for the repulsive 
term, leading to very different estimates, depending on the chosen model (see the table). 
Therefore it is tedious to predict well depths with an accuracy better than a few tens of wave 
numbers. An empirical solution to this problem is offered when an experimental determina- 
tion of the ground state potential is available, as illustrated for example in the spectroscopic 
analysis of NaRb states by Docenko et al |28j|: the difference between the experimentally 
determined ground state potential curve and the computed curves of ref.J3] is attributed 
to the cumulative effect of the different terms of the core-core interaction. This difference is 
then assumed to be independent of the molecular electronic excitation: it can be added in 
turn to the computed excited potential curves, yielding then a good representation of the 
excited molecular states. 
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Atom 


basis (A) 


basis (B) 


basis(C) 


Li 


8s7p4dlf/[7s5p3dlf] 


10s7p5d3f 




Na 


7s5p5d2f/[5s5p3d2f] 


8s6p5d2f 




K 


7s5p7d2f /[6s4pAd2f] 


7s5p7d2f 




Rb 


7sAp5dlf /[6s4p4dl/] 


9s6p6d4f 




Cs 


7s4p5dlf/[Qs4p4:dlf} 


7s4p5d4f 


9s6p6d4f 



TABLE I: Gaussian basis sets for each studied alkali atom. The contractions are specified in 
brackets for the "A" basis. The "C" basis is set up only for Cs. 
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TABLE II: Exponents of the Gaussian functions introduced 
in the various basis A, B, C. When appropriate, contracted 
orbitals are displayed in brackets, with contraction coeffi- 
cients in parenthesis. 



Atom 


Basis 


i 


Exponents 


Li 


A 


s 

P 
d 
f 

s 

S 

f 


2.464(-0.0f3829),f.991(-0.032077)J,0.582,0.2,0.07,0.03f, 0.015,0.007 
0.630(0.052433),0.240(0.078566)],0.098,0.043,0.02,[0.01(0.01445),0.005(0.00657)] 
0.2,0.07,[0.022796(0.491521),0.008574(0.549272)] 
0.01 

2.464,1.991,1. ,0.582,0.2,0.136,0.07,0.011,0.008,0.004 

0.630,0.24,0.098,0.03,0.02,0.01,0.004 

0.2,0.062,0.024,0.01,0.004 

0.01,0.0125,0.04 




B 


Na 


A 
B 


s 

P 

d 
f 
s 
P 

A 

f 


2.8357(0.007043),0.49318(-0.1871),0.072085(0.3679)j,0.036061,0.016674,0.00693,0.00287 
0.431(-0.01778),0.09276(0.2003)],0.03562,0.01447,0.0058,0.00023 
0.292(0.01454),0.06361(0. 123)], 0.02273, 0.008852,0.00352 
0.015,0.0055 

2.8357,1. ,0.49318,0.072085,0.036061,0.017674,0.00693,0.00287 

0.431,0.09276,0.03562,0.1447,0.0058,0.00023 

999 06361 09973 008859 00359 

0.015,0.0055 


K 


A 
B 


s 

P 

d 
f 

s 

P 
d 
f 


0.9312(0.02463),0.2676(-0.2628)|,0.0417,0.02815,0.01448,0.0055,0.0026 
0.133,0.05128,0.01642,0.0052,0.0022 

1.255(0.02754),0.4432(0.05391),0.109(0.1083)],0.02994,0.01013,0.0037,0.0018 
0.015,.005 

0.9312,0.2676,0.051,0.02715,0.01448,0.0055,0.0026 
0.134,0.05328,0.01642,0.0053,0.002 
1.255,0.443,0.108,0.02994,0.01013,0.0037,0.0018 
0.015,0.005 


Rb 


A 


s 

P 
d 

f 

s 

P 

d 
f 


1.292561(-0.065450),0.824992(0.037862)|,0.234682,0.032072,0.013962,0.005750,0.0025 
0.128,0.040097,0.014261,0.004850 

0.408807(0.0943),0.096036(0.1846)],0.026807,0. 009551,0. 0004 




B 


2:5,1-292561,0.824992,0.5,0.234682,0.032072,0.013962,0.00575,0.025 
0.25,0.128,0.040097,0.03,d.014261,d.00485 
0.328807,0.2,0.096036,0.026807,0.009551,0.0004 
0.2,0.1,0.05,0.005 


Cs 


A 


s 

P 
d 

f 

s 

!! 

f 

s 

P 

d 
f 


0.328926(0.411589),0.241529(-0.682422)|,0.050502,0.029302,0.013282,0.00528,0.003 
0.12,0.0655,0.0162,0.00443 

0.196894(0.18965),0.067471(0.22724)],0.027948,0.010712,0.003 

8:328926,0.221529,0.050502,0.029302,0.013282,0.00528,0.003 
0.155,0.0655,0.0162,0.00443 
0.18894,0.067471,0.027948,0.010712,0.003 
0.05,0.025,0.0125,0.005 

2.5,0.328941,0.241529,0.1,0.050502,0.029302,0.012,0.005,0.002 

0.25,0.12,0.0655,0.03,0.012,0.005 

0.3,0.196894,0.09,0.03 

0.05,0.025,0.0125,0.005 




B 

C 
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Atom 


Basis 






Pp 


Pd 


Pf 


Li 


A 


0.19151451 


1.434 


0.979 


0.6 


0.4 




B 


0.1997[73j 


1.315 


0.999 


0.6 


0.9 


Na 


A 


0.9947[45l 


1.42 


1.625 


1.5 


1.5 




B 


0.9987[73j 


1.456 


1.467 


1.5 


1.5 


K 


A 


5.354[45l 


2.067 


1.905 


1.96 


1.96 




B 


5.472[73j 


2.11439 


1.97355 


1.9884 


1.9884 


Rb 


A 


9.245 


2.513 


2.279 


2.511 


2.511 




B 


9.245[£3j 


2.5538 


2.349 


2.5098 


2.5098 


Cs 


A 


15.117[2U 


2.6915 


1.8505 


2.807 


2.807 




B 


16.33^,2^ 


2.8478 


1.981 


2.904 


2.904 




C 


16.33[73, 75] 


2.8864 


2.7129 


2.8963 


2.8963 



TABLE III: Dipole polarizabilities and semiempirical cutoff parameters pi introduced in the 
basis A, B, C. For previously published basis (A), the authors used the theoretical values quoted 
in the quoted references for a% 1+ , while for basis (B) we used instead the experimental values of 
Wilson et al Note that for Cs, Allouche and coworkers 
do not actually use his value. 
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76J referred to Wilson et al, but 
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Aa 


A s 






Aa 


A s 






Aa 


As 


Li 


2s 


5.67 


0.55 


Na 


3s 


-2.26 


0.00 


K 


4s 


0.10 


-0.18 




2p 


0.21 


-0.76 




3p 


-1.46 


0.47 




4p 


0.03 


0.14 




3s 


7.79 


26.67 




4s 


19.81 


2.50 




5s 


-5.49 


-9.99 




3p 


1.80 


8.47 




3d 


0.64 


0.01 




3d 


0.58 


-0.61 




3d 


7.94 


0.20 




4p 


3.94 


2.67 




5p 


47.01 


39.38 




4s 


1684.04 


54.70 




5s 


11.75 


7.54 




4d 


23.07 


21.30 




4p 


4069.42 


236.05 




4d 


13.70 


13.39 




6s 


3.01 


0.01 




4d 




40.36 




4f 




57.45 




4f 


46.32 


46.30 




4f 


721.84 


51.28 




6p 


50.12 


59.90 














A A 


A s 






A A 


As 


Ac 




Rb 


5s 


0.08 


-0.03 


Cs 


6s 


0.32 


0.15 


0.31 










5p 


0.00 


0.00 




6p 


0.05 


-0.85 


0.10 










4d 


0.13 


0.03 




5d 


-18.85 


-0.51 


0.46 










6s 


-32.17 


-32.56 




7s 


1.50 


10.40 


7.09 










6p 


32.67 


22.92 




7 P 


68.43 


91.63 


17.40 










5d 


7.29 


4.15 




6d 


186.58 


203.32 


85.30 










7s 


-10.39 


-10.79 




8s 
4f 


37.55 


41.02 
23.64 


20.20 
23.47 









TABLE IV: Differences (in cm 1 ) between computed and experimental atomic binding energies 
for A, B, and C basis sets. Experimental values are taken from refs. 
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Molecule 


LiNa 


LiK 


LiRb 


LiCs 


NaK 


NaRb 


NaCs 


KRb 


KCs 


RbCs 


basis (A) 


98 


104 


89 


89 


114 


99 


99 


105 


105 


90 


basis (B) 


142 


148 


162 


149 


136 


143 


137 


156 


143 


157 


basis (C) 








162 






156 




156 


170 



TABLE V: Size of the molecular basis set generated from method A, B and C. 
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TABLE VI: Computed permanent dipole moments \D\ of the 
X l Y* + state (in Debye) at equilibrium distances R£ , at the 
minimum distance R$ (in a.u.), and for the v = level of the 
X state. Our results are compared to available experimental 
and theoretical values. 





Tl 


lis work 


Exp. 


Th 




D{Rf) Rf 


D(Rf) Rf 

1 Li LL 


D(v = 0) 


D R e ref. 


D(R e ) R e ref. 


LiNa 


(A, 


561 5 43 
0.554 5.42 


0.630 6.85 
0.633 6.89 


0.566 
0.556 


0.45( ) 162 
0.47(3) 5.33 
0.463(10) 5.31 M 
0.45 ' f65 


1.24 5.64 156 
0.485 5.429 
0.485 5.42 g7 
0.487 5.47 g£ 

0.53 5.42 [52' 


LiK 


(A, 
(B) 


3.558 6.21 
3.533 6.21 


3.807 7.50 
3.792 7.49 


3.565 
3.555 


3.510(5) 
3.45(2) 6.18 
3.87 


63 
61 
65 


3.437 6.292 |45j 
3.50 6.25 


LiRb 




K 

B 


4.168 6.52 
4.142 6.48 


4.442 7.78 
4.414 7.78 


4.165 
4.131 


4.05 


65j 


4.13 6.52 [52] 


LiOs 




A 
B 

X! 


5.520 6.81 
5.512 6.82 
5.462 6.81 


6.023 8.33 
5.998 8.31 
5.970 8.33 


5.529 
5.524 
5.478 


6.30 [65J 


5.48 6.89 [52] 


NaK 




A 


2.760 6.50 
2.763 6.49 


2.854 7.42 
2.862 7.47 


2.759 
2.762 


2.73(9) 
2.690(14) 6.55 
3.42 


64 
61 
65 


3.6 6.9 
2.735 6.36 
2.95 6.614 
2.54 6.59 
2.78 6.6 


54 
45 
55] 
52] 
35 


NaRb 




:| 


3.304 6.84 
3.301 6.84 


3.413 7.75 
3.413 7.76 


3.306 
3.301 


3.10(3) 6.73 
3.51 


§1 
65 


3.33 6.85 


52] 


NaCs 




A 
B 

'& 


4.613 7.20 
4.661 7.20 
4.580 7.20 


4.821 8.26 
4.864 8.26 
4.793 8.29 


4.607 
4.660 
4.579 


4.75(20) 6.91 
4.75 
5.86 


61 
78' 
65 


4.60 7.23 [52] 


KRb 




A 

;b) 


0.615 7.64 
0.589 7.64 


0.620 8.20 
0.605 8.49 


0.615 
0.589 


0.20 


65j 


0.64 7.65 152 
0.71 7.7 M 
1.06 7.7 [48' 


KCs 




A, 

B 

C 


1.906 8.01 
1.921 8.02 
1.835 8.02 


1.957 8.87 
1.967 8.85 
1.891 8.93 


1.906 
1.921 
1.837 


2.58 [65] 


1.92 8.05 152] 


RbCs 




A 


1.238 8.28 
1.278 8.30 
1.205 8.30 


1.278 9.19 
1.309 9.12 
1.240 9.19 


1.237 
1.280 
1.204 


2.39 [65] 


1.26 8.71 [52] 
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TT-disp 
Vcc 


[3JJ 


[S21-A [S2J-E 


KRb 


7.64 


-6.41 




«17 ~-507 


RbCs 


8.30 


-10.22 


18.42 


«25 «-500 



TABLE VII: Contribution (in cm 1 ) of the Vcc' sp and Vlc P terms around the equilibrium distance 
of the KRb and RbCs ground state. The Vcc. sp term is estimated follo wing equation ^ and 



parameters from method A. Both two results (labelled as A and E) from ref. 69] are displayed. 
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V. FIGURE CAPTIONS 



Figure H {Color on line.) The potential curves V'(R) (including only the 1/R term 
of V CC (R)) around their minimum for (a) the and (b) the a 3 S + states of RbCs, as 

obtained through the three approaches A (black full lines), B (red dashed lines), C (blue 
dashed line with dots), described in the text. In (a), results for B and C approaches are 
superimposed. In addition, the black dot-dashed lines represent V'(R) + V™ d (R) for both 
states, in the A case. 

Figure El (Color on line.) The permanent electric dipole moments (in Debye) of (a) the 
X 1 S + (b) the a 3 E + states, in RbCs, as obtained through the three approaches A, B, C, 
described in the text. 

Figure H2 (Color on line.) The permanent electric dipole moments (in Debye) of (a) the 
X 1 S + state (b) the a 3 £ + state, as functions of the internuclear distance, obtained from 
calculation B. Results are displayed in black for Li coumpounds (dashed line for LiNa, 
dot-dashed line for LiK, double-dot-dashed line for LiRb, full line fo LiCs), in red for Na 
compounds (open circles for NaK, upper triangles for NaRb, closed circles for NaCs), and in 
blue for K compounds (lower open triangles for KRb, crosses for KCs) and for RbCs (lower 
closed triangles). 

Figure EJ (Color on line.) The permanent electric dipole moments (in Debye) of the 
X X S+ state of the mixed alkali pairs as functions of the vibrational level, obtained from 
calculations B. The color and line code is the same than in figure El 

Figure El (Color on line.) The permanent electric dipole moments (in Debye) of (a) the 
A 1 S + state (b) the a 3 S + state of KRb as functions of the internuclear distance, obtained 
from calculations A (full line) and B (dashed line). Our results are compared to other 

and relativistic (closed 
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4j| (dot-dahed blue 



available theoretical determinations: non-relativistic (open circles) 
circles) results from ref.j3], and non-relativistic results from ref. 
line). The vertical dashed line in (b) panel indicates the location of the repulsive wall of the 
potential at the energy of its dissociation limit. 
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Figure 1, Aymar and Dulieu, Journal of Chemical Physics 
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Figure 2, Aymar and Dulieu, Journal of Chemical Physics 
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Figure 3, Aymar and Dulieu, Journal of Chemical Physics 
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Figure 4, Aymar and Dulieu, Journal of Chemical Physics 
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Figure 5, Ay mar and Dulieu, Journal of Chemical Physics 




